Effective Temperature Characterizing the Folding of a Protein 
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The time sequences of the molecular dynamics simulation for the folding process of a protein is 
analyzed with the inherent structure landscape which focuses on configurational dynamics of the 
system. Time dependent energy and entropy for inherent structures are introduced and from these 
quantities an effective temperature is defined. The effective temperature follows the time evolution 
of a slow relaxation process and reaches the bath temperature when the system is equilibrated. We 
show that the nonequilibrium system is described by two temperatures, one for fast vibration and 
the other for slow configurational relaxation, while the equilibrium system is by one temperature. 
The proposed formalism is applicable widely for the systems with many metastable states. 

PACS numbers: 05.70.Ln 87.15. Cc 05.90.+m 



Slow relaxations of non-equilibrium systems are a very 
general process which are typically observed in glasses. 
The focus of the study for glasses has been put on descrip- 
tions of their instantaneous state, or their aging proper- 
ties 0, often characterized by an effective temperature 
deduced from the violation of the fluctuation-dissipation 
theorem On another hand, for various slow re- 

laxations a major interest is how the relaxation proceeds. 
This is the point that we want to investigate in this let- 
ter by showing how an effective temperature that follows 
the time evolution of a slowly relaxing non- equilibrium 
system can be introduced. The effective temperature in 
glasses is, in most cases, defined from the ratio between 
equilibrium and non-equilibrium response to aperturba- 
tion, and therefore it is observable-dependent [g . On the 
contrary, the effective temperature that we introduce in 
this letter is a one-time property of the system, which 
evolves towards the equilibrium temperature as the sys- 
tem relaxes. 

To illustrate the usefulness of the concept, we have 
chosen to consider proteins, which fold from random coil 
states to their own native conformation in /is or ms, while 
their vibrational dynamics takes place in time scales of 
the order of /s or ps. The folding of a protein provides a 
nice example of an out-of-cquilibrium complex physical 
system which evolves in an experimentally and numeri- 
cally accessible time scale. Proteins can be investigated 
in the spirit of glasses by characterizing their inherent 
structures, i.e. the local minima of their energy landscape 
• We showed that the density of states in the Inherent 
Structure Landscape (ISL) of a model protein depends 
exponentially on their energy 0, so that the model has 
some similarities with the simple models of glasses such 
as spin glasses 8], Lennard- Jones clusters 9] and the 
trap model . This letter shows that the ISL is not only 
useful to study the equilibrium properties of a protein T*] 
but can also be used in practice to characterize its slow 
evolution towards the folded configuration. 

Model: We consider an off-lattice Go protein model 
with a slight frustration . The protein is reduced to its 



backbone without side chains, each residue being repre- 
sented by its Ca atom. The effective potential between 
the Ca elements in the backbone is designed from the 
experimentally determined native structure. We study 
proteinG which has 56 residues 0|. At low temperature, 
the model protein folds to the native structure, while at 
high temperature it denatures. Folding-unfolding transi- 
tion occurs at Tf, and around Td ~ 0.4Ty a dynamical 
transition is observed below which the fluctuations of the 
protein structure are dramatically reduced. 

In numerical investigation, molecular dynamics simu- 
lations with Langevin equations are applied, where a de- 
natured initial condition is prepared by a simulation at 
1.7 Tf at which the protein is in the random coil state. 
A sudden cooling to T <Tf is applied at i = 0. For each 
value of T, 200 initial conditions are sampled. 

Definitions of quantities: Slow dynamics comes from 
itinerant transitions among local energy well while the 
fast dynamics comes from rapid vibrations inside each 
well. Such slow dynamics can be characterized by "in- 
herent structures" j in which the set of conformations 
belonging to the same basin labeled by the index a are 
mapped to its local energy minimum with energy Cq. 

Let ^(r, t, T) be the probability density for the system 
to be at point F, up to dT, in phase space at time t and 
let Ta be the region corresponding to a basin a. The 
probability weight to observe the basin a is 



Wait, T) 



w{T,t,T)dT, 



(1) 



where X)a^"(^'-^) ^ ^- "^^g probability weight 
w'^{T ,t,T) normalized in the basin a is introduced as 



,{T,t,T)^Wa{t,T)-w':{V,t,T). 



(2) 



To be consistent with Eq.©, /^^ u;^(r, t, T)dT = 1. 

For each point F, the value of the energy e(r) is de- 
termined from the model potential. When the point 
F belongs to basin a, the energy can be rewritten as 
e(F) ^ Ca + AVq(F), where AVq(F) is the potential at 
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point r measured from the minimum Cq. The mean en- 
ergy U at time t is given by 



U{t,T) EE J e{T)w{T,t,T)dr 

= Uis{t,T) + J2Kit,T)wa{t,T). (3) 

This expression splits the energy into two compo- 
nents, Uis{t,T) = ^„e„«;„(t,T) and U^it,T) = 
/pA14(rX(r,t,r)dr. Uis{t,T) is the mean inherent 
structure energy at time t which reflects slow structural 
changes due to the transition among the basins while 
U"{t,T) is the mean vibrational energy inside basin a. 

Since the inherent structure energies Cq, take discrete 
values, the entropy for the inherent structures can be 
defined with the definition of quantum statistical physics, 

Sis {t,T) = -kBj2^o.(t, T) log {t,T). (4) 

a 

The total entropy S{t,T) splits into Sis{t,T) and a vi- 
brational contribution S"{t,T) according to 

s{t,T) = -fee y w(r,t,r)iogu;(r,i,T)dr 

= 5is(t,r)+^5,"(t,r)u;„(t,T). (5) 
Here S"{t,T) is the vibrational entropy at the basin a 

such that s^{t,T) = -kB /p^ <(r,t,r)iog<(r,t,r). 

We add that so-called configuration entropy is not nec- 
essarily equal to 5*18 even in canonical equilibrium. 



^IS,can(r) ^ fcBfogas(t/lS,ca„(r)), 



(6) 



because the value of Bq, can largely deviate from the mean 
value Uis^ca,n{T) in small systems such as proteins. Here 
^^13(60) is the density of states for the inherent struc- 
tures, which can be determined from a numerical simu- 
lation in equilibrium situation as shown in 0- 

In analogy to the equilibrium case, we define a time- 
dependent effective temperature for inherent structures 
as 
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_ dSis{t,T)/dT 

TAt,T) - dUis{t,T)/dT' 



(7) 



where we expect Te(t, T)^T as the relaxation proceeds, 
and Te = T in thermal equilibrium. 

Slow relaxation in Inherent structure landscape: The 
energetic aspects of the folding process are demon- 
strated in Fig^ in which the mean vibrational en- 
ergy over the phase space is defined as Uv{t,T) = 
J2aUyit,T)wa(t,T). As are shown in Fig^ the mean 
inherent structure energy Uis{t,T) relaxes logarithmi- 
cally in time. Such logarithmic relaxations are typically 
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FIG. 1: Time evolution of mean inherent structure energy Uis 
and (inset) mean vibrational energy Uv In the case 0.55T/ > 
Td, the folding process reaches thermal equilibrium around 
t — 2 X 10^, whereas in O.SATf < Td it does not within the 
computation time. Uv is in the equilibrium value which is 
slightly larger than the harmonic approximation i/2k^NrT 
(indicated as 0.51 and 0.83 for T/Tf = 0.34 and 0.55). 



seen in glassy systems, which suggests intrinsic glassy 
properties for proteins even above the dynamical transi- 
tion temperature T^j. 

On the other hand, U^{t,T) is almost constant, cor- 
responding to its equilibrium value at temperature T, 
although the probability weight Wa {t, T) is expected to 
significantly depend on time as the folding proceeds. This 
is possible if the vibrational energy U^{t,T) is approxi- 
mately the same for all basins a and the vibrational mo- 
tions reach equilibrium much faster than the time scale 
of the folding. The value of S^{t, T) is also expected to 
be equal to its equilibrium value, which could be approx- 
imately independent of a, similarly to U^(t,T). 

Sig{t,T) cannot be accessible to numerical simulation 
because the number of basins in the whole phase space 
is too large to determine Wa{t, T) for each a and each t. 
Instead of determining Wa{t,T) for all basins, let us in- 
troduce a probability P(eQ,i,T) to observe the inherent 
structure energy Cq, 

P(ea, t,T) = J dTwir, t, T)(5(e,„i„(r) - e„). (8) 

Here einin(r) is the quenched energy from the point T. If 
r belongs to the basin a, ei„in(r) — e^. 

Since the initial conditions in numerical simulation are 
randomly selected from the random coil state, let us 
assume that the weight is approximately the same for 
basins with the same minimum energy Cq, 



y,(t,T) ^P(e„,t, r)Ms(e„). 



(9) 



Then, 



Sisit, T) ^ -fcB /de„P(e„, i, T) log (10) 
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FIG. 2; Time evolution of the probability distribution 
P{ea,t,T) for the inherent structure. The horizontal axis 
is inherent structure energy Ca scaled by ksTf. T = O.SATf 
and 0.557/ correspond to T < Td and T > Td, respectively. 
The right bottom figure shows the equilibrium distribution at 
T = 0.55T/. At t = the protein has been cooled to T. 
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FIG. 4: (a) Time evolution of Te{t,(Ti + T2)/2) in 
eq.(|7J estimated by {Sis{t,T2) - Sis{t,Ti)} / {T2 - Ti) and 
{Uis[t,T2) - Uis{t,Ti)}/{T2 ~ Ti). Here (Ti.Ta) = 
(0.347/, 0.4ir/),(0.4ir/, 0.487/) and (0.557/, 0.597/). The 
dotted horizontal line is (7i + 72)/2 for (0.557/, 0.597/), 
which is expected to be a convergence line of 7e(t, 0.577/). 
(b) Confirmation of thermodynamic relation in Eq. (??) in 
thermal equilibrium condition. 
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FIG. 3: Time evolution in the plot on Sis vs Uis for t > 10*". 
7/7/ = 0.34 and 0.55. The line is the relation between Sis 
and Uis in thermal equilibrium. Higher temperatures corre- 
sponds to the right-upper part of the line. 
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FIG. 5: (a) lJ^{t,T) vs Uis{t,T). Here T = {Ti + 
72)/2 with (7i,72) = (0.347/, 0.417/), (0.417/, 0.487/) and 
(0.557/, 0.597/). 



which is accessible to numerical simulations. 

Examples of the time evolution of P{ea,t, T) in Fig[21 
show that the distribution is strongly dependent on time 
after the sudden cooling. The energy range that the dis- 
tribution spans becomes smaller as time proceeds, and 
the width of the distribution decreases significantly. As 
is expected P(eQ,,i,T) is considerably different from the 
equihbrium distribution (bottom figure for T/T/ — 0.55). 

For each time t, we can obtain numerically Uis(t,T) 
and Sis{t,T) following the definition I|1U|) . If we plot 
Sis{t,T) as a function of Uisit^T) for various t, we get 
the time evolution for the relaxation process shown in 
FigOl For comparison we plot the line corresponding to 
thermal equilibrium cases over a wide range of T. This 
comparison makes it clear that the time evolution in a 
fixed cooling temperature T proceeds in parallel to the 
line for thermal equilibrium from high to low tempera- 



tures. Let us study this evolution by the effective tem- 
perature Te{t,T) defined in eq. Q. In equilibrium. Sis 
and Uis work well to determine the temperature of the 
system as is demonstrated in Fig. 0Jb), where Tg for the 
canonical ensemble with ^is{ea) has the expected prop- 
erty, Te = T. In the relaxation process, Tg decreases 
with time from a higher temperature. When T is suffi- 
ciently high, we can observe that Te converges to a value 
corresponding to T after the relaxation to equilibrium 
(T ~ 0.55T/ in Fig.l^a)). Thus the features observed in 
Fig. m or Fig. O in the relaxation process are reflected in 

Te(i,T). 

The value of Te(i, T) is expected to give an approxima- 
tion of P{ea,t,T) to an equilibrium distribution at the 
temperature Tg. To check it, we compare the approxi- 
mated mean IS-energy Ue{t,T) with Uisit,T) in Fig.|Sl 
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where 

(11) 

The figure displays the tendency Ue{t,T) ~ Uis{t,T), 
which suggests that Tg can give a first approximation for 
the distribution P{ea,t,T). The deviation of the points 
from the fine in Fig. O indicates that the states on the 
way of folding do not coincide strictly with any equilib- 
rium state. Correspondingly, the points in Fig.|21deviates 
slightly from the equilibrium line. 

Discussion: We have introduced time dependent en- 
ergies and entropies for inherent structures to analyze the 
slow folding process of a model protein. After a sudden 
cooling, we found that the probability distribution of the 
inherent structures is widely spread in the early stage 
and becomes narrower as time proceeds. This evolution 
is successfully characterized by the time evolution of the 
effective temperature Eq.(j7|l, as shown in Fig0fa). 

Although the approximation of the distribution by 
is not complete, we note that even this first approxima- 
tion is not obtained if we treat the total probability distri- 
bution without separating the inherent structure modes 
and the vibrational modes. As argued in Fig^ the vibra- 
tional modes are reaching local equilibrium in the early 
stage. Thus there is a temperature difference between 
the fast dynamics and the slow dynamics. Thermody- 
namic frameworks including two temperatures T and 
could be proposed [ij]. It might be interesting to study 
various slow dynamics adopting such two temperatures. 

When T is lower than Td, relaxation becomes very slow 
as is expected, but it would relax after a long period 
because the present model is small finite system without 
extensibity. We expect a convergence tendency of to 
T with time. This could be in contrast to spin glasses, 
since the shape of P{ea, t, T) for a spin glass [3 does not 
change in the slow relaxation process. Then Te might 
be obtained as a constant larger than T, which could be 
compared with the effective temperature in the violation 
of fluctuation-dissipation system. 

The effective temperature defined in the present pa- 
per is different from the known effective temperature in 
glasses. For instance, Sciortino and Targlia introduced 
an effective temperature according to the inherent struc- 
ture energy and configuration entropy in Lennard- Jones 
systems. They successfully showed that their effective 
temperature coincides with the one defined by the viola- 
tion of fluctuation-dissipation relation. These effective 
temperatures are the parameter to indicate the prop- 
erty of the individual system. Even though the system 
is very close to equilibrium, their effective temperature 
does not converge to the surrounding temperature. On 
the other hand, the effective temperature of Eq.(|2j) is de- 
fined to indicate the relaxation process of the system. 



It becomes equal to the surrounding temperature if the 
system reaches equilibrium. 

Finally we would like to point out that this study em- 
phasizes the interest of the Inherent Structure Landscape 
for protein. We showed earlier that, contrary to the 
free-energy landscape which is only a useful concept, the 
density of states of the ISL can be obtained in practice 
for a protein model, and that it characterizes the equilib- 
rium properties of the protein to a good accuracy. Here 
we show that the ISL is not only useful in equilibrium, 
but can be used to characterize the slow relaxation with 
a time-dependent effective temperature. The assumption 
used in this letter, that and 5*" does not depend on 
a, appears to be good for the protein case. It can be 
relaxed for more complex systems jl3|. 

Note: Since the present model possesses at least two 
funnel which is well separated from the others, we choose 
the relaxation process to a funnel of which global mini- 
mum energy state corresponds to a misfolded state. Its 
minimum energy state resembles the native state but a 
ck-helix is bent compared to the native state. In this 
funnel, relaxation proceeds in a time scale accessible to 
numerical simulations if the temperature is not too low. 
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